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ABSTRACT 

Using a sub-pixel event repositioning technique, we spatially resolved X-ray emission from the 
infrared double system IRS 5 in the R Corona Australis molecular cloud with ~0."8 separation. 
As far as we know, this result - obtained from 8 Chandra archival observations between 2000 
and 2005 - is the first X-ray study of individual sources in a Class I protostar binary system 
with a projected separation of less than 200 AU. We extracted light curves and spectra of the 
individual sources using a two-dimensional image fitting method. IRS 5a at the south, the source 
which was brighter in the near-infrared, showed three X-ray flares lasting >20 ksec, reminiscent 
of X-ray flares from pre-main sequence stars, while the northern source (IRS 5b) was quiescent 
in X-rays in all the observations except for a 2005 August 9 observation with a factor of ^2 
flux enhancement. In quiescence, these sources showed almost identical X-ray spectra, with JVh 
~4xl0 22 cm~ 2 , &T~2 kcV, and log Lx ^30.2—3 erg s _1 . IRS 5a showed plasma at temperatures 
up to kT ~5— 6 kcV during flares, while the column density of IRS 5b increased by a factor of 
2 during an observation on 2005 August 9. We discuss the evolutionary stages and variation of 
the X-ray activity of these sources. 



Subject headings: techniques: high angular resolution 
stars: pre-main-sequence — X-rays: stars 

1. Introduction 

Many, if not most, stars are born as a mem- 
ber of a multiple system through fragmentation of 
the parent molecular cloud. Because young stellar 
objects (YSOs) that form as a binary or multiple 
system are of the same age and chemical compo- 
sition, any differences in the level of their stellar 
activity should be due to differences of stellar ini- 



binaries: visual — stars: magnetic fields 



tial parameters such as mass, net angular momen- 
tum and/or the presence of the companion star 
(See reviews about multiplicity of youn g stars in 



Duchene et alj|2007c iMonin et al.ll2007l and mag 



netic activity on various stellar parameters in ex. 



Chabrier fc Kukerll2006l ; iDobler et al.ll2006l ). YSO 
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binaries are therefore a useful probe to test models 
of the dependence of the activity on stellar param- 
eters for stars of the same age. 

The R Corona Australis (R CrA) molecular 
cloud has rich association of YSOs with vari- 



ous evolutionary stages (jTavlor fc Storevl 119841 ). 
Among them, IRS 5 (a.k.a. TS 2.4) suffers strong 
extinction (Ay ^ 45 mag) and is clas s ified as a 
Class I protostar (jWilking et all Il986l . Il992i see 
therein for the spectral classification of YSOs). 
High angular resolution near-infrared imaging sep- 
arated IRS 5 into two stellar components, IRS 5a 
(K ~10.9 m ) and IRS 5b (K ~11.5 m ). The sep- 
aration of the two stars is less than an arcsec- 



ond (|Chen fc Grahamlll993l : lNisini et al.ll2005l ). In 
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what follows the term IRS 5 means the entire bi- 
nary system, or is used when the emitting source 
was unknown due to limited imaging capability. 
The southern source, IRS 5a, has a stellar spec- 
tral type of K5-K7V and is located i n a low mass 
proto star phase on the HR diagram ( Nisini et al.l 
20051 ). Stellar parameters of IRS 5b are poorly 
known due to its faintness. 

IRS 5 is also known as one of the rare Class 
I protostars showing st rongly variable cm con- 
tinuum radio emission (ISuters et all 1 19961) with 



signifi cant circular polarization (jFeigelson et al 
1998). Such gyrosynchrotron radio emission is 



believed to be the major radio emission mecha- 
nism of weak-lined T-Tauri Stars (wTTSs, e.g. 
Carkner et al.lll997l ). but the emission from pro- 
tostars is supposed to be absorbed by the sur- 
rounding ionized materials heated by mass accre- 
tion. IRS 5 does not currently show strong mass 
accretion activity, but there is eviden ce of an out 



flow p ossibly driven by IRS5 (HH 73 1. 1 Wang et al 
2004). The near-infrared spectrum of IRS 5a 



shows small continuum excess above the photo- 
spher ic flux, which also suggests a low accretion 
rate (jNisini et al.l 120051 ) . The infrared spectrum 



also shows charact eristics of scattere d emission. 
From these results. iNisini et al. I (|2005l) questioned 



the idea that the IRS 5 system is a bona-fide pro- 
tostar, but it may be close to the TTS phase. 
IRS 5 also showed strongly absorbed thermal 



Kovama et al. 


19961: Neuhauser & Preibisch 19971 


Forbrich et al. 


20061. 20071 Forbrich & Preibisch! 



panied by occasional rapid flares. These X-ray 
characteristics are consistent with plasma heating 
by magnetic activity in the stellar corona, so that 
the X-ray and radio emission from IRS 5 is sus- 
pected to_arise_jrom the sa me root cause. Sim- 
ilarly, iForbrich et al. ( 2007 ) did not find a clear 



correlation between the X-ray and near-infrared 
flux variation from IRS 5 and other cluster mem- 
bers, which tends to support a coronal origin for 
the X-ray activity instead of accretion driven ac- 
tivity. 

To understand the relation between the X-ray 
and radio activity of IRS 5 and study the evolu- 
tion of the activity, we need to resolv e the IRS 5 
binary system in both wavelengths. I Choi et af 



radio sources in a VLA radio observation. Re- 
solving X-ray sources with a sub-arcsecond sepa- 
ration is challenging even with the Chandra ob- 
servatory with its unprecedented X-ray angular 
resolution. However, encouraged by a marginal 
resoluti on of IRS 5 in a Cha ndra image (see Fig- 
ure 5 of IForbrich et alJ 120071 ). we revisited Chan- 
dra archival observations of IRS 5 and clearly 
resolved the individual binary components using 
enhanced spatial information obtained with the 
Sub-pixel Even t Repositioning (SER) tec hnique 
(|Tsunemi et al.1 120011 iLi et al.l EM 12004 ). Us- 
ing this technique we were also able to extract 
light curves and spectra of each source with a two- 
dimensional image fitting method and measured 
the plasma parameters of individual binary stars 
for the first time. 

The distance to the c loud is controversial. 
Marraco fc Rvdgre"r] ( 1981 ) estimated it to be 
129 pc from uvby/3 photometry of t hree stars in 
the cl oud. Assuming Ay =3.1 mag. ICasev et al 



(|1998I ) derived the same distance (129±11 pc) to 
the eclipsing binary TY CrA, whi ch is suspected 
to be a cluster member. However, iKnude fc H(bq 
(jl998l ) measured d = 170 pc from the interstellar 
reddening of stars in the cloud having HIPPAR- 
COS distances and Tycho B - V colors, and whose 
distance we use for IRS 5. 

2. Observations 

We found eight arch ival pointings with Chandra 
(jWeisskopf et al.ll2002t ) that include IRS 5 in the 
FOV (Table [l}. These observations occurred be- 
tween 2000 October and 2005 August and used the 
Advanced CCD Imaging Spectrometer with the 
Imaging array (ACIS-I) at the prime focus with ex- 
posure times between 15—40 ksec. During the ob- 
servations, IRS 5 was put at 1'— 2' off-axis, where 
the PSF quality is negligibly different from the 
on-axis position. We use the convention that indi- 
vidual Chandra observations are designated CXO, 
subscripted with the year, month and day of the 
observation. We downloaded the data from the 
HEASARC archive^ and reprocessed them with 
the CIAC@ analysis software version 3.4 (CALDB 
ver. 3.3.0.1). We reprocessed the data, removed 



(|2008l ) successfully resolved IRS 5 into two cm 



1 http: / /heasarc. gsfc.nasa.gov/W3Browse/ 
2 http:/ /cxc. harvard.edu/ciao/ 
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Table 1 

Journal of the Chandra Observations 



Image Shift 



Abbreviation 


Seq. ID 


Observation Start 


Exposure 


A (R.A., Dec.) 


StdDev 


Off-axis 








(ksec) 


(arcscc) 


(mas) 


(arcscc) 


CXO001007 


19 


2000 Oct 7 17:02 


19.7 


(0.40, 0.96) 


128 


59.3 


CXO030626 


3499 


2003 Jun 26 12:58 


37.6 


(0.56, 0.28) 


81 


90.0 


CXO040617 


4475 


2004 Jun 17 23:17 


19.9 


(0.19, 0.21) 


117 


114.2 


CXO 50808 


5402 


2005 Aug 8 2:38 


15.2 


(0.07, 0.06) 


88 


72.1 


CXOososos 


5403 


2005 Aug 9 2:39 


15.0 


(0.06, 0.11) 


109 


72.4 


CXO050810 


5404 


2005 Aug 10 1:58 


15.0 


(0.08, 0.06) 


67 


72.5 


CXO050812 


5405 


2005 Aug 12 3:13 


15.0 


(0.06, 0.12) 


144 


72.5 


CXO050813 


5406 


2005 Aug 13 1:51 


14.6 


(0.03, 0.01) 


18!) 


72.6 



Note. — Seq. ID: sequence identification number of each observation. Image Shift: shift of the X-ray 
sky coordinate frame from the original X-ray data to the 2 MASS frame. StdDev: standard offset 
deviation of the X-ray source positions from the 2MASS source positions in milli-arcsecond. Off-axis: 
off-axis angle of IRS 5 from the nominal point. 



the pixel randomization and utilized the SER al- 
gorithms to improve the spatial resolution. 

Photon pile-up can cause source count rates 
to be underestimated, spectral shapes to appear 
harder, and degrades event position determination 
in the SER method. Fortunately, observed X-ray 
photon counts of IRS 5 ranged between 0.02—0.08 
counts s , corresponding to small pile-up frac- 
tions of ^2—8 9(Q We thus did not correct for 
any pile-up effect in our analysis. 

3. Analysis and Results 
3.1. Absolute Astrometry 

This study crucially requires absolute astrome- 
try of the X-ray images with a few tenths of arc- 
second precision. We therefore measured sky po- 
sitions of X-ray sources around IRS 5 and cross- 
correlated them with sources in the Point Source 
Catalog of the 2MASS All-Sky Data Release. 

There is up to 7 year time span between the in- 
terval when the 2MASS observations of the south- 
ern sky (1998 March - 2001 February) and the 
Chandra observations (2000 Oct - 2005 Aug). 
The cluster members can move significantly within 
7 years if their proper motions are suitably large. 
The HIPPARCOS catalogue h as 4 sources with 
prop er motion measurements (jPerrvman fc ESA 
19971 ). Among them, R CrA and V709 CrA have 
proper motion of ~80— 100 mas yr _1 , but they 



have large uncertainty in the parallax measure- 
ment (~68, 140 mas) and almost opposite direc- 
tions in their proper motions. This would suggest 
that their measurements arc less reliable, possi- 
bly due to contamination of surrounding nebu- 
lar emission. Two other sources, TY CrA and 
HD 176386, with better parallax measurements 
have small proper motions (<20 mas yr _1 ). Our 
tentative measurement of the proper motion of 
cluster members using VLA data between 1997 
and 2005 also derived s imilarly small pr oper mo- 
tions of <28 mas vr^ ljChoi et al.l l2008h . These 
results suggest that the proper motion of the clus- 
ter members is small, less than ~150 mas in 7 
years (<1 pixel in the Figure 1 image). 

To measure X-ray source positions in each ob- 
servation, we ran the wavdetect source detection 
tool in the CIAO package for a 0.35—8 keV image 
binned at 0."5 pixel -1 and only considered sources 
at above the 4ct significance. The 2MASS posi- 
tions of IRS 1, IRS 2, R CrA, IRS 13 and HBC 667 
had corresponding X-ray counterparts with simi- 
lar offsets (<0."1 bet ween the sources, see sou rce 
positions in Figure 1 of lHamaguchi et al. I l2005l and 
Forbrich et al.ll2006l ). while there was no counter- 



part to the X-ray source corresponding to IRS 6 
in the 2MASS catalogue within 1"0 We there- 
fore used the former 5 sources for the positional 
references. 

The 2MASS positions of IRS 1, IRS 2, R CrA, 



3 http: / /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



4 It was in the "reject" table. 
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Fig. 1. — Magnified images of IRS 5 (raw: left, smoothed: right), sequentially in time from above. The gray 
scales are shown in a linear scale for all figures (but saturated above 10 cnts for the left panels). Crosses 
show X-ray source positions determined from our analysis. Circles at the bottom left of the left panels show 
the half of the peak intensity of the PSF. Raw image bins are 0."123 pixel -1 . Gray scales and contours in 
the right panels are normalized by the effective area and exposure time in CXOq.sosos- 
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Fig. 1. — Continued. 
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Fig. 2. — Composite of images in CX0 oioo7) 
CXO 05 0808, CXO 05 08i2 and CXO 0508 i 3 , when both 
of IRS 5a and IRS 5b were in the quiescent 
phas e. Markers — cr oss: X-rays, circle: cm ra- 
dio (|Choi et all |2008[) . diam ond: 2MAS S , star : 
infrared source positions in iNisini et all (|2005l ). 



shifted to match for the cm radio source posi- 
tions. Error ranges of the X-ray source positions 
from the 2MASS fram e (200 mas), radio beam 
size (IChoi et al. 2008!) and IR image pixel size 
(jNisini et al.ll2005h arc shown at the bottom left. 



IRS 13 and HBC 667 (K -7.0-10.5" 1 ) have pos- 
sible position inaccuracy of up to 100 mas RMS, 
except for the brightest source R CrA (K =2.9 m ), 
which had a slightly larger position uncertainty 
of 160 mas (RMSjE The uncertainty of the mean 
position of all the sources is ~50 mas. Com- 
bined with the cluster proper motion noted above, 
we can determine the absolute source position 
better than ^200 mas for these sources. Ta- 
ble [T] lists the shifts of the X-ray sky coordi- 
nate frame from the original X-ray data. Ear- 
lier observations in 2000—2003 needed relatively 
large shifts (~0.3— 1"), while observations in 2005 
needed almost negligible shifts (<0.2"). Table Q] 
also lists the standard offset deviation of the X- 
ray source positions from the 2MASS source posi- 
tions. The result suggests that X-ray source posi- 
tions in all the Chandra observations were deter- 



mined at <200 mas from the 2MASS frame, and 
hence <400 mas in the absolute coordinate frame. 

3.2. Image Analysis 

Figure Q] shows X-ray images of IRS 5 between 
1—8 keV after the absolute position correction. 
They were binned with 0."123 pixel -1 to resolve 
sources having sub-arcsecond separation. In all 
the images, X-ray emission was clearly elongated 
toward the NE-SW direction. We checked in two 
ways if this structure is real. First, we created 
simulations of point spread fun ctions (PSF) at the 
location of IRS 5 with ChaRT (jCarter et al.ll2003l ) 
and MAR]^ version 4.2.1, assuming a spectrum iVn 
=4.5xl0 22 cm -2 , kT =3 keV and Z =0.3 solar, 
(based on our preliminary fits of the IRS 5 spectra 
using the standard method). Second, we made im- 
ages of IRS 1 with a similar off-axis angle to IRS 5 
to determine if the SER method might introduce 
asymmetries in the PSF. We found no such exten- 
sion from either of the half maximum PSF circles 
(Figure Q]) nor the IRS 1 images. We conclude 
that the deduced elongation of the X-ray image of 
IRS 5 is real. 

In CXOo30626j CXO040617 and CXOososio the 
south-west part of the elongation was remark- 
ably bright, while in CXO050809 the north-east 
part looked slightly brighter than the south-west 
part. In the other observations, we recognized 
apparently two peaks with similar intensities. We 
re-ran the wavdetect source detection tool for 
these high resolution images and detected the 
south-west source in all observations, and detected 
the north-east source in CXOooioo7j CXOo50808> 
CXO050809 and CXOo508i2- The north-east peak 
was not detected when the south-west peak was 
very bright or images lacked sufficient expo- 
sure. Both peaks were located at, on average, 
(a 20 oo, £2000) = (19 /l l m 48. s 07, -36°57'21."8) 
designated using the Chandra naming conven- 
tion as CXOU J190148. 1-365722 and (a 20 oo, 
£2000) = (19 h l m 48. s 02, -36°57'22."4) as CXOU 
J190148.0-365722 and the positions of the indi- 
vidual detections differed by 0."09 for the former 
and 0."075 for the latter in RMS. These sources 
are separated by only ^0."8±0."12 (ler), corre- 
sponding to a ^140 AU projected separation at 
d = 170 pc. 



3 http: / / www. ipac.caltech.edu / 2mass / releases /allsky/doc/ 



3 http://space.mit.edu/CXC/MARX/ 



6 



QO 



FO 



FO 



QO 



FO 



FO 



QO 



QO 



'( ( 

+ IRS5a 




If 




( ( 1 ' ( 


( 1 ' ( ( 


1 ' ( 


( 1 1 




- +IRS5b 






+ 






+ 
+ 










+ 


+ 




+ 


+ 

+ 










Hz 




















H — h 










001007 


030626 


it 

040617 


it — H 

050808 


r 1 — 

050809 


050810 


r — — / 

050812 


050813 



Observations 



Fig. 3. — Light curves between 1—8 keV (top) of IRS 5a (black) and IRS 5b (red). Each bin has 4 kscc and 
tic marks on the horizontal axis are shown by 10 ksec from the beginning of each observation. The error bars 
show ler confidence range. In CXO030626, CXO040617, CXO050810, and CXO050813, IRS 5b emitted X-rays 
comparable to the other observations, but wavdetect did not detect it. This is probably because enhanced 
X-ray emission from IRS 5a contaminated strongly at the IRS 5b position, and thus increased the apparent 
background level (see M3.3[) . Labels at the top depict observations which showed significant X-ray variations 
in either of IRS 5a nor IRS 5b (FO: Flare Observation) or in neither of them (QO: Quiescent Observation). 



A composite image of CXO001007, CXO0508O8, 
CXO050812 and CXO050813 clearly showed two 
peaks at the north-east and south-west sides 
(Figure \2§ ■ These X-ray peaks have correspond- 
i ng cm radio cou nterparts within ~0."1 for each 
(|Choi et al.ll2008l ). They have offsets from the 
infrared sources IRS 5a and IRS 5b of ac- 
cording to their absolute coordinat es in the H 



and K band maps in Figure 1 of iNisini et al 
(|2005h . The 2MASS survey with an absolute 
astrometric accuracy of ~100 mas detected an 
infrared source very close to the north-east X-ray 
source though it did not resolve the two infrared 
sources probably due to the limited spatial reso- 
lution of the CCD (2"). Moreover, the two X-ray 



(200E 



sources and the IR sources in INisini et al. 
have similar position angle and separation, and 
therefore match their positions well by shifting 
the IR image frame (see Figure [2]). We therefore 
identified the south-west X-ray source (CXOU 
J190148.0-365722) as IRS 5a and the north-east 
source (CXOU J190148.1-365722) as IRS 5b. 



3.3. Timing and Spectral Analysis with 2- 
Dimensional Image Fits 

Because the X-ray emission from IRS 5a and 
IRS 5b is strongly blended, the standard method 



to generate light curves and spectra — count- 
ing photons within user defined source and back- 
ground regions — cannot be used straightfor- 
wardly. Therefore we utilized the 2-dimensional 
image fitting method in which we generated multi- 
ple images from each observation in restricted time 
and energy ranges and fit each 2-dimcnsionally by 
the combined individual model PSF images for 
IRS 5a and IRS 5b. We created images with a size 
of 80 pixel x 80 pixel binned by 0."123 pixel -1 . 
For the model images, we used the ptsrc2d func- 
tion in the CIAO package, assuming two point 
source components at the position of the X-ray 
sources, and used the sherpa tool to carry out 
the fitting procedure (We should note that a part 
of the X-ray emission may come from different 
and/or extended regions.) We did not include 
background in the model images, since the back- 
ground is neglig ibly small (~ 3xl0~ 4 cnts / (1 
pixcb=0."123) 2 / 10 kscc). The ptsrc2d func- 
tion loads an image as a template model PSF. 
Using these model PSFs, we generated photon 
events at the posit ion of IRS 5 using both ChaRT 
([Carter et al.ll2003l ) and MARX. For the light curve 
analysis, we assumed a PSF to have the typ- 
ical spectrum of IRS 5 based on our prelimi- 
nary analysis using the standard method, i.e. iVn 
=4.5xl0 22 cm -2 , fcT=3 keV and 0.3 solar elemcn- 
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Fig. 4. — Quiescent spectra of IRS 5a (top left) and IRS 5b (top ri<?/i£) and flare spectra of IRS 5a in 
CXO030626 (middle left), CXO040617 (middle right), CXOososos (bottom left) and of IRS 5b in CXO050809 
(bottom right), overlaid on the quiescent spectra (grey). 
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Table 2 
Spectral Fits 



Source 


Spectrum 


Photon 


Model 


kT 




X 2 /d.o.f (d.o.f.) 








(counts) 




(keV) 


(10 22 cm" 2 ) 




(10 30 ergs" 1 ) 


IRS5a 


QOa 


490.6 


IT 


2.2 (1.7-3.2) 


3.7 (3.0-4.5) 


0.78 (24) 


1.5 




FOao30626 


771.9 


IT (+QOa) 


5.4 (>2.8) 


4.7 (3.3-6.8) 


0.58 (64) 


2.5 




FOao406i7 


554.1 


IT (+QOa) 


5.4 (>2.8) 


5.1 (3.7-6.8) 


0.67 (52) 


4.4 




FOaososio 


597.5 


IT (+QOa) 


6.6 (>3.4) 


5.1 (4.1-6.5) 


0.69 (57) 


6.2 


IRS5b 


QOb 


493.0 


IT 


2.0 (1.5-2.5) 


4.5 (3.7-5.8) 


0.92 (24) 


1.9 




FObo50809 


292.0 


IT (+QOb) 


4.4 (>1.4) 


7.8 (5.1-14.4) 


0.89 (43) 


3.3 



Note. — QO: Quiescent Observation. FO: Flare Observation. Photon: net photon counts between 1—8 keV. Absorption 
corrected Lx between 0.5 — 10 keV, assuming the distance of 170 pc. Values in parentheses denote formal 90% confidence 
limits, and do not include any systematic errors produced by our 2-dimcnsional fitting procedure. The elemental abundance 
was fixed at 0.3 solar value. 



tal abundance. In both the timing and spectral 
analyses, we generated PSF images in the same 
energy band as the observed images using large 
numbers of photon events to minimize statistical 
fluctuations. In the sherpa fits, we fixed the cen- 
troids of two point source components at the X-ray 
source positions of IRS 5a and IRS 5b and only al- 
lowed their normalization parameters to vary. We 
used the Cash statistic appropriate for Poisson- 
limited data and used the Powell method to find 
the best-fit values. We derived upper and lower 
boundaries of the \a confidence range using the 
projection function. 

Figure [3] shows individual light curves of 
IRS 5a and IRS 5b in 4 ksec bins in the en- 
ergy range 1—8 keV. We normalized the individ- 
ual lightcurves to the 3 keV effective area at the 
position of IRS 5a, resulting in small correction 
factors (<1 %). Their flux variations look to be 
uncorrelated, suggesting no significant interfer- 
ence with each other in our 2-dimensional fits. To 
further check the consistency of the image fitting 
results, we summed the individual light curves of 
IRS 5a and IRS 5b and compared them to a light 
curve extracted from the entire IRS 5 complex, 
including photon counts within a 3" radius cir- 
cle centered on between IRS5a and IRS5b, which 
include ~90— 93% of photons from both IRS 5a 
and IRS 5b. We found that the standard method 
obtained 93—98% of the photon counts we derived 
by combining the individual IRS 5a and IRS 5b 
photons. We had similar results for the split en- 
ergy bands between 1—3.5 keV and 3.5—8 keV, 
which we used for a hardness ratio analysis de- 



scribed in the next paragraph. This suggests that 
our image fitting method produces a result that is 
consistent with the standard method within <10% 
discrepancy. 

In CXOo30626i CXOo406i7i and CXO050810, 
IRS 5a showed significant time variations, reminis- 
cent of stellar X-ray flares. The fast rise (~4 ksec 
in CXO040617) and slow decay (a half decay time 
of 15—30 ksec) time scales based on visual in- 
spection were similar to those of relatively strong 
X-ray flares from YSOs. The flux level did not 
change remarkably during the other observations 
(<50%). IRS 5b was quite stable in almost all 
the available observations. The non-detection of 
IRS 5b with wavdetect in 4 observations would be 
therefore caused by increased contamination from 
IRS 5a by flares, or from limited photon statis- 
tics. IRS 5b showed a small flux enhancement 
by a factor of 2 in CXO050809, and an apparent 
decrease at the middle of CXO040617, but this 
decline could be an artifact produced by strong 
contamination from the IRS 5a flare in this obser- 
vation along with a possible offset of the south- 
west emission peak from IRS 5a. We split the 
energy band into two at 3.5 keV, which is near 
the median of the photon distribution against 
energy in a spectrum, and measured the hard- 
ness ratio (HR) of each observation, which we 
have defined as HR = (H - S)/(H + S), where 
the H and S are count rates of the 3.5—8 keV 
and 1—3.5 keV bands, respectively. The average 
HRs of both IRS 5a and IRS 5b in CXO001007, 
CXO0508O8, CXO050812 and CXO050813, when nei- 
ther showed significant variations (and which we 
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refer to as the "quiescent observations"), were 
HR ~ -0.32±0.15, -0.29±0.18 (la), respec- 
tively. The HRs were 0.02±0.09 (la) dur- 
ing flares from IRS 5a (CXO 30626, CXO 406i7, 
CX0 05 o8io) and ~0.04±0.08 (la) during a small 
flux enhancement of IRS 5b (CXOo50809)- The 
HR was slightly higher in the flare than during 
the quiescent observations. 

We generated individual spectra of IRS 5a and 
IRS 5b in key phases and analyzed them with 
xspec version 11. One technical issue was that 
xspec requires integer photon counts, while the 
ptsrc2d function outputs the normalization val- 
ues (i.e. photon counts) in floating point. Thus, 
to avoid any round off error, we multiplied pho- 
ton counts and exposure time by 1000 when cre- 
ating spectral files. The xspec format also as- 
sumes a single error value for each spectral bin, 
while the projection function outputs upper and 
lower error ranges. We therefore used the larger 
value output from projection as the xspec error 
for each spectral bin. Since we set most of the 
spectral bins to have >15 photons the upper and 
lower error ranges typically have less than ~10% 
differences. We generated response matrices and 
auxiliary files corresponding to the spectral files 
using mkacisrmf and mkarf in the CIAO package. 

Since no quiescent observation had enough pho- 
ton statistics to generate a spectrum with suffi- 
cient signal-to-noise, we combined images of these 
observations in each spectral bin. We combined 
their PSF images and spectral responses as well for 
our 2-dimcnsional image fits and spectral model 
fits. The result is shown in Figure 21 The qui- 
escent spectra of IRS 5a and IRS 5b were almost 
identical with significant emission up to ~8 keV 
and a cut-off due to absorption at ~1.5 keV. We 
modeled both spectra using an absorbed optically 
thin, 1-temperature (IT) plas ma model (WABS, 
Morrison fc McCammonlll983t APE(f|, Table E). 
The best-fit models had Nn ~4xl0 22 cm -2 and 
kT ^2 keV for both sources. The absorption cor- 
rected log L x was 30.2-3 (0.5-10 keV). 

We compared the IRS 5a spectra during flaring 
states (CXO 30626, CXO 406i7 and CX0 05 o8io) to 
its quiescent spectrum (Figure^). All flare spectra 
showed stronger emission in the 1—8 keV energy 
band than the quiescent spectrum. We fit the flare 



spectra by an absorbed IT plasma model plus the 
best-fit model of the quiescent IRS 5a spectrum, 
assuming that a new plasma component appeared 
during the flares. The best fit to the flare spectra 
seems to have higher kT (~5— 7 keV) than the qui- 
escent spectrum but similar JVh (~5x10 22 cm -2 ). 
However, IRS 5b showed similar increase in hard 
band flux in CXO050809 compared to the quies- 
cent phase, while the flux below ^2 keV was ap- 
parently unchanged (Figure 2]). The best fit, as- 
suming an absorbed IT model plus the best-fit 
model of the quiescent IRS 5b spectrum showed 
significantly larger N K -7.8l*; 2 xl0 22 cm" 2 (la), 
which does not overlap with iVjj in the quiescent 
spectrum (4.5±^ x 10 22 cm~ 2 , la). The best-fit 
model also showed slightly higher kT ~4 keV than 
the quiescent spectrum. None of these spectra had 
enough photon statistics to show the Fe line pro- 
file, which appeared in some flare spectra of the 
entire IRS 5a complex. 

The values of iV H and fcTfor IRS 5a and IRS 5b 
are similar in all the observations and a re con - 

(|200fih . 



sistent with the values in Forbrich et al 



which were derived from the combined spectra of 
IRS 5a and IRS 5b. The kT and Lx values are 
typical of those from TTSs and Class I protostars, 
while the derived Nyj is typica l of Class I proto- 



stars (e.g.. Ilmanishi et al.ll200lh . 



4. Discussion 



We still have a discrepancy betwe en the posi- 
tions of the near-infrared sources in iNisini et al 



(|2005l ) and the sources detected by 2MASS, VLA 
and Chandra. IRS 5a was a fa c tor of two brighter 
than IRS 5b in INisini et all (|2005f) . while the 



2MASS source, whose position should be biased 
toward the brighter source (IRS 5a), was instead 
closer to IRS 5b (see fl3~2|) . The 2MASS source is 
unlikely to be mis-positioned due to surrounding 
contamination since both IRS 5a and IRS 5b are 
su fficiently bright in t he images shown in Figure 1 
of INisini et aTT (|2005h . IRS 5a or IRS 5b might 



7 http:/ /cxc. harvard.edu/atomdb/ 



have experienced an increase in its near-infrared 

(2005) 



flux during the observations in INisini et al. 
or during the 2MASS observations. We should 
note that IRS 5 only showed variation of less than 
±0.03mag duri ng monitoring observa tions in 2005 
(see Figure 6 of iForbrich et al. 2007 ). 

It is possible that the X-ray sources we have 
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Table 3 

Young Binary System with <200 AU Projected Separation Resolved in X-rays 



Object 


System 


Separation 


Distance 


Age 


Reference 






sky 


projection 












(arcscc) 


(AU) 


(PC) 


(Myr) 




TWA 5 


M3Ve (TTS) + M8.5-9 (BD) 


2 


110 


55 


12 


Tsuboi et al. (2003) 


HD 100453 


A9Ve + M4-M4.5 (TTS) 


1.06 


114 


120 


20 


Collins et al. in prep. 


HD 98800 


K5Ve (TTS) + ? 


0.8 


48 


38 


5-10 


Kastner et al. (2004) 


R CrA IRS 5 


K5-K7V (Class I) + ? 


0.8 


140 


170 


0.3-0.5 or 1.3 a 


this result 



Note.— TTS: T-Tauri Star, BD: Brown Dwarf. 
"Values of IRS 5a from lNisini et al.l 1 12003) . 



identified with either IRS 5a or IRS 5b could be 
a background red giant star or an active galac- 
tic nucleus, though this is unlikely. However, X- 
ray emission from red giant stars are not usu- 
ally as bright as, nor as hot as, emis sion from 
IRS 5a or IRS 5b (|Pizzolato et al.ll2000h . Though 



some AGNs are observed with fluxes similar to 
those from IRS 5a or IRS 5b, the X-ray spec- 
tra are ge nerally harder (s pectral photon index 
T <2, ex. 



Ueda et al. 1998h than those observed 
for either IRS 5a or IRS 5b (T ~3.5 assuming an 
absorbed power-law mod el). On the other hand, 



Pontoppidan et all (|2003h showed CO ice absorp- 
tion lines of IRS 5a and IRS 5b have similar sys- 
temic velocities to the other YSOs in the R CrA 
cloud (~ —20 km s _1 ) though this feature can be 
produced by an absorber physically unrelated to 
the infrared sources. The X-ray spectra are typi- 
cal X-ray emission from Class I protostars. These 
results suggest that IRS 5a and IRS 5b form a 
protostar binary system within the R CrA cloud. 

Table [3] shows binary systems with projected 
separations less than 200AU that have been spa- 
tially resolved in X-rays. According to this table, 
IRS 5 is the first X-ray binary system believed to 
be comprised of Class I protostara^ and is a fac- 
tor of >10 younger than the other systems. Bi- 
nary systems with linear separations of less than 
140 AU have very little disk mass, perhaps because 
of the g ravitational influence of the binary com- 
panion (jBeckwith et al.lll990h . The IRS 5 binary 
separation is near this boundary and the separa- 



3 Although the infrared spectral slope, which defines the 
YSO classification, has not been measured for IRS 5b, we 
think this star is a Class I protostar based on its high ab- 
sorption in X-rays. 



tion may be smaller if the distance is d ~130 pc. 
The separation is so large that the binary stars 
would not have direc t interaction by magne tic 
field, for example (see Uchida fc Sakurail 19851 as 
an example of the system with magnetic interac- 
tio n), but they may have an influence as discussed 
by iBeckwith et al. ( 1990h . IRS 5 may be an im- 
portant example of the influence on X-ray activity 
of binary interaction in a very early phase. 

The column densities of IRS 5a and IRS 5b 
derived from their X-ray spectra (~ 4x 10 22 cm" 2 ) 



were typical of Class I protostars (c.f. llmanishi et al 
200 ll ). These values of Nn were a factor of 2 
smaller than the JVh value estimated from Ay 
for IRS 5a (~45 mag), assuming an empirical 
Nh— Ay relation app ropriate to the p Oph cloud 



( Imanishi et al. 200ll ). This result was consistent 
with X-ray absorption me asurements of IRS1, 2 
and 5 ( Forbrich et al. 20061 ) and qualitatively con- 



sistent witli_a_sjnallJ^r/Aj ratio of the R CrA 
cloud (|Vuong et al.ll2003fh The moderate N H may 



not suggest the e dge-on geometry discussed by 



Nisini et alj (|2005D since X-rays from the stellar 



core would suffer stronger extinction in the edge- 
on disk than infrared emission from the stellar 
core and inner disk. 

The kT and Lx values of the quiescent spectra 
of both IRS 5a and I RS 5b are typical of C lass I 
protostars and TTSs (jlmanishi et al.l 120011 ). As- 
suming the mass o f IRS 5a between 0.4—0.9 M Q 



(jNisini et al.ll2005D . the L x ~10 30 2 erg s _1 corre 



sponds to a ~1 M-year-old pre-main sequence star 
based on the Lx-age relation of young stars in the 
Orio n nebula though the relation has a large scat- 
ter (jPreibisch fc Feigelsonl [iooih . The kTs dur- 
ing the flares of IRS 5a are typical of YSOs (e.g., 
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Imanishi et~ai1l200ll ; IWolk et alj|2005l) . 

IRS 5a exhibited 3 prominent flares during the 
8 Chandra observations. Since the observations 
ended at the middle of these flares, we only know 
lower limits to the flare durations (>37.6, 19.9 and 
15.0 ksec) and total flare energies (>9.2, 8.7 and 
9.3 xlO 34 ergs), respectively. These flares were 
similar in brightness to that of the brightest flares 
from classical TTSs and weak -line TTSs observe d 
in the XEST project (Fig. 9 of lStelzer et ail2007t) . 
Curiously, we did not see any flare-like event from 
IRS 5a fainter than those flares though such flares 
would have been detectable. IRS 5a had 3 flares 
during the 152 ksec exposure, that is, one flare per 
~50 ksec. This is more than an order of magni- 
tude more frequent than the average of the XEST 
(1 flare in 770 kse c for flares with the flare energies 
above 10 35 ergs, IStelzer et al] 12007ft and COUP 
(1 flare in 650 ksec. IWolk et al .1 120051 ). though we 
did not consider possible detections of flares that 
start outside the observing window. IRS 5 showed 
remarkable variat ion in the other XM M-Newton 
observations too ( Forbrich et al. 20061 ). for which 



either of IRS 5a or IRS 5b would be responsi- 
ble. IRS 5a seems to be an actively flaring source, 
though the accreti on rate was measured to be low 
(jNisini et al.ll2005l ). 

IRS 5b did not show apparent variation dur- 
ing the Chandra observations except for a flux en- 
hancement in CXOo50809- The column density of 
IRS 5b during the flux enhancement increased by 
a factor of ~2 from the quiescent phase. Such 
variable X-ray emission in CXO050809 might have 
been produced from a plasma on or around the 
stellar surface that appeared behind the circum- 
stcllar disk, or perhaps IRS 5b may have another 
deeply embedded companion. 

Except for the flare activity, X-ray properties 
of IRS 5a and IRS 5b arc strikingly similar. This 
may strengthen a general assumption that stars 
in a similar mass range experience similar evolu- 
tion of their quiescent X-ray spectra. This result 
also raises a caution in measurements of the Lx 
function of YSOs — Lx can be significantly over- 
estimated if a YSO has an unresolved companion. 

5. Conclusion 

We spatially resolved X-ray emission from each 
component of the infrared double system IRS 5 in 



the R CrA molecular cloud with a 0."8 projected 
separation in 8 Chandra observations. Their X-ray 
brightnesses and spectra suggest that they are a 
Class I protostar binary system. Separations of bi- 
nary YSOs are typically less than ~10" in nearby 
star forming regions. Chandra enabled study of 
X-ray properties of such young binary systems for 
the first time. 

We derived light curves and spectra of these ob- 
jects using a 2-dimensional image fitting method. 
During the 8 Chandra observations, the southern 
source IRS 5a showed three X-ray flares lasting 
more than 15 ksec and with flare energies of >10 35 
ergs. These flares were relatively powerful com- 
pared to flares from classical TTSs and weak-line 
TTSs. The flare frequency (1 flare per ~50 ksec) 
was higher than flares from PMSs in the Taurus 
and Orion cloud. IRS 5a might be in an active 
phase, or may have a different mechanism which 
drives its X-ray flares. IRS 5b, on the contrary, 
did not show strong variation through the obser- 
vations, except for a flux increase by a factor of 
two in 2007 Aug. 9 accompanied by an apparent 
TVh increase. 

Though IRS 5a and IRS 5b had clear differences 
in their X-ray variability, their quiescent X-ray 
spectra were almost identical. This may suggest 
that both stars have gone through a similar evolu- 
tionary phase in X-ray activity. Detailed infrared 
spectroscopy of IRS 5b is crucial to determine its 
stellar parameters and to understand nature of the 
IRS 5 system. Measurement of proper and radial 
motions of both IRS 5a and IRS 5b would be also 
important to understand whether they arc a grav- 
itationally bound system. 

This work is performed while K.H. was sup- 
ported by the NASA Astrobiology Program un- 
der CAN 03-OSS-02. This publication makes use 
of data products from the Two Micron All Sky 
Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and 
Space Administration and the National Science 
Foundation. This research has made use of data 
obtained from the High Energy Astrophysics Sci- 
ence Archive Research Center (HEASARC), pro- 
vided by NASA's Goddard Space Flight Center. 

Facilities: CXO(ACIS-I) 
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